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Abstract. Low ice crystal concentration and sustained in- 
cloud supersaturation, commonly found in cloud observa- 
tions at low temperature, challenge our understanding of cir- 
rus formation. Heterogeneous freezing from effloresced am- 
monium sulfate, glassy aerosol, dust and black carbon are 
proposed to cause these phenomena; this requires low up- 
drafts for cirrus characteristics to agree with observations 
and is at odds with the gravity wave spectrum in the upper 
troposphere. Background temperature fluctuations however 
can establish a “dynamical equilibrium” between ice produc- 
tion and sedimentation loss (as opposed to ice crystal forma- 
tion during the first stages of cloud evolution and subsequent 
slow cloud decay) that explains low temperature cirrus prop- 
erties. This newly-discovered state is favored at low tempera- 
tures and does not require heterogeneous nucleation to occur 
(the presence of ice nuclei can however facilitate its onset). 
Our understanding of cirrus clouds and their role in anthro- 
pogenic climate change is reshaped, as the type of dynamical 
forcing will set these clouds in one of two “preferred” micro- 
physical regimes with very different susceptibility to aerosol. 


1 Introduction 

Cirrus clouds are composed of ice crystals that form at 
high altitudes and temperatures typically below 235 K (Prup- 
pacher and Klett, 1997). They play a key role in climate 
by modulating the planetary radiative balance (Liou, 1986) 
and heat transport in the upper troposphere (Ramanathan and 
Collins, 1991). They strongly impact water vapor transport 
across the tropopause level (Jensen and Pfister, 2004) and 
play an important role in lower stratospheric chemistry (Pe- 
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ter, 1997). Cirrus may be affected by aircraft emissions (Se- 
infeld, 1998) and long range transport of pollutants and play 
an important (but highly uncertain) role in anthropogenic cli- 
mate change. 

A key microphysical parameter required for understand- 
ing the climate impact of cirrus is their concentration of ice 
crystals, N c . At temperatures between 200 and 235 K cirrus 
ice crystals form primarily by homogenous freezing of super- 
cooled deliquesced aerosol (DeMott et al., 2003; Heymsfield 
and Sabin, 1989), which occurs if the saturation ratio with 
respect to ice, .S’, (i.e., the ratio of water vapor partial pres- 
sure to its equilibrium value over ice) reaches a characteris- 
tic threshold value. Shorn (Koop et al., 2000). Heterogeneous 
freezing of water upon existing aerosol particles (termed “ice 
nuclei”, IN) can also occur (at S lower than Shorn) an d con- 
tribute to ice crystal concentrations (DeMott et al., 2003; 
Froyd et al., 2009), especially in polluted and dust-rich re- 
gions (Barahona et al., 2010a; Haag et al., 2003). The level 
of water vapor supersaturation (i.e., S — 1) is the thermody- 
namic driver for ice formation, and is generated by expansion 
of air parcels forced by large scale dynamics, gravity waves, 
and small scale turbulence (Kim et al., 2003). 

At temperatures below 200 K (typically near the tropical 
tropopause layer, TTL) the simple conceptual model for cir- 
rus formation presented above is at odds with observations 
(Jensen et al., 2010; Kramer et al., 2009; Peter et al., 2006). 
Temperature fluctuations from mesoscale gravity waves are 
common at high altitudes and can produce localized vertical 
motion with updraft velocity as large as 1 ms -1 (Bacmeis- 
ter et al., 1999; Herzog and Vial, 2001; Jensen and Pfister, 
2004; Sato, 1990). Such motion can increase the rate of ex- 
pansion cooling at the point of freezing so that a large num- 
ber of ice crystals is nucleated before the local supersatura- 
tion is depleted by ice crystal growth. Homogeneous freez- 
ing driven by gravity wave motion would produce high ice 
crystal number concentration, N c , between 1 and 10 cm -3 
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near the TTL (Barahona et al., 2010a). Such high concentra- 
tions however are not observed; A c remains low, sometimes 
even lower (0.005-0.2 cm -3 ) than concentrations observed 
in weak updraft zones at cold temperatures (Kramer et al., 
2009; Lawson et al., 2008). This “low A c ” paradox is ac- 
companied by other phenomena, such as low supersaturation 
relaxation times (Kramer et al., 2009), which in turn leads to 
sustained supersaturation levels inside clouds (i.e., “the su- 
persaturation puzzle”) (Gao et al., 2004; Kramer et al., 2009; 
Peter et al., 2006), high clear-sky supersaturation (Jensen et 
al., 2005), and broad ice crystal size distributions (i.e., large 
crystal sizes, Jensen et al., 2008). These phenomena occur 
despite the strong dynamical forcing and the ample amounts 
of deliquesced aerosol available for homogeneous freezing. 
Suppressed freezing by organics (Murray, 2008), slow wa- 
ter vapor transfer to the ice phase (Gao et al., 2004; Magee 
et al., 2006), and freezing to cubic instead of hexagonal ice 
(Murray et al., 2005), have been proposed to explain these 
features. These mechanisms are however not capable of ex- 
plaining the low A c in low temperature cirrus clouds (Peter 
et al., 2006). Lacking the predictive understanding of such 
phenomena hinders the ability of climate models to capture 
the climate effects of cirrus clouds and their response to an- 
thropogenic perturbations. 

Heterogeneous freezing of IN as the main path of cirrus 
formation has been proposed to explain the features of cir- 
rus clouds at low temperature (Abbatt et al., 2006; Jensen 
et al., 2010; Murray et al., 2010). Owing to their ability to 
freeze ice at much lower supersaturation than homogeneous 
freezing requires, IN can deplete water vapor, reduce super- 
saturation and inhibit homogeneous freezing; this can drasti- 
cally reduce the number of ice crystals that forms in the cirrus 
(Barahona and Nenes, 2009b; DeMott et al., 1994; Karcher 
et al., 2006). Much of the anthropogenic impact on cirrus 
clouds and climate is thought to occur through this IN-A C 
feedback mechanism (Lohmann and Feichter, 2005). Dust 
(Khvorostyanov et al., 2006), effloresced ammonium sulfate 
(Abbatt et al., 2006; Jensen et al., 2010; Wise et al., 2010), 
and glassy aerosol (Murray et al., 2010) have been identified 
as potential heterogeneous IN at the TTL . 

The evolution of cirrus clouds at low T has been addressed 
in several studies. Using a one dimensional (1-D) cloud 
model, Jensen and Pfister (2004) found that the superposi- 
tion of temperature fluctuations along Lagrangian trajecto- 
ries near the TTL resulted in rapid cooling cycles that in- 
creased the rate of crystal production by homogeneous nu- 
cleation leading to A c above 1cm -3 . Khvorostyanov et 
al. (2006) used a 1-D cirrus model to investigate the evo- 
lution of a cirrus layer initialized at 200 K and found that 
regardless of the predominant nucleation mechanism, tem- 
perature fluctuations increased the maximum A c in the cloud 
(up to 0.6 cm -3 ). A c however rapidly decreased after the ini- 
tial freezing pulse due to the vertical advection of ice crys- 
tals precluding new nucleation events and diluting A c down 
to about 0.05 cm -3 . Gensch et al. (2008) used box model 


simulations along Lagrangian trajectories to test homoge- 
neous and heterogeneous nucleation scenarios in the forma- 
tion of cirrus at low T . It was found that only heterogeneous 
nucleation scenarios (with prescribed IN number concentra- 
tion around 0.1 cm -3 ) would result in A c close to observa- 
tions. This conclusion was echoed by Froyd et al. (2009) and 
Jensen et al. (2010) using 1-D models along prescribed T 
trajectories. 

The formation of cirrus clouds however exhibits complex 
non-linear behavior that may not be captured by box and 1 - 
D models. Ice falling through active freezing zones (typi- 
cally located at the top of the cirrus layer, Spichtinger and 
Gierens, 2009b) in clouds consume water vapor and can in- 
hibit homogeneous freezing much like IN do (Kay et al., 
2007; Spichtinger and Gierens, 2009b). Their effectiveness 
depends on their residence time in freezing zones and hence 
depends on their size. Large ice crystals tend to quickly fall 
out of freezing zones and have limited effect on new ice for- 
mation events; small crystals (typically those with terminal 
velocity, M term , less or equal to the mean updraft u of the 
cirrus layer) fall slowly and can remain long enough in the 
upper part of the cloud to affect new freezing events. This 
suggests that at low temperatures, preexisting (and typically 
small, Kramer et al., 2009) ice crystals may locally dehydrate 
the freezing zone sufficiently to inhibit the formation of new 
ice. The rate of crystal production is not uniform through the 
freezing zone, as the “local” saturation ratio, .S', and updraft 
velocity, u (defined at the scale of individual cloud “parcels” 
~10°-10 2 m, Pruppacher and Klett, 1997) may be affected 
by fluctuations in wind speed and temperature induced by 
gravity waves (Jensen et al., 2010; Karcher and Haag, 2004; 
Kim et al., 2003). These internal S variations are usually ne- 
glected in cirrus cloud studies on the basis that the long-term 
evolution of the cloud is determined by the mean values of S 
and u. 

In this work we analyze the range of conditions for which 
heterogeneous freezing may explain the features of cirrus 
clouds at low temperature, and propose an alternative view 
(based on a statistical description of cirrus formation and 
evolution) in which the interplay of temperature fluctuations, 
and ice crystal production and sedimentation leads to pre- 
viously unidentified natural cirrus states of low ice crystal 
concentration and sustained high supersaturation. 

2 Heterogeneous freezing at low temperature 

The impact of IN on A c depends on their concentration, 
Ain- If too low (Ain < 1 x 10 -4 cm -3 ), a negligible im- 
pact is seen on A c , as too few (heterogeneously-frozen) 
ice crystals form to quench supersaturation below the ho- 
mogeneous/heterogeneous freezing threshold (Barahona and 
Nenes, 2009b). Low A c favors large crystal size and there- 
fore heterogeneously frozen ice crystals may sediment out of 
the cloud layer before significantly modifying S (Spichtinger 
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and Gierens, 2009a). When Ain approaches a characteristic 
“limiting” concentration (which depends on updraft velocity, 
the IN freezing threshold and size), Au m , supersaturation is 
quenched, homogeneous freezing is suppressed, and A c de- 
creases steeply (Barahona and Nenes, 2009b). For Ain > 
A]i m , homogeneous nucleation is inhibited and A c = Ain- 
Thus, Au m is the minimum A c that can form in an active 
nucleation zone in a freshly-formed cirrus cloud (Barahona 
and Nenes, 2009b) and presents the maximum reduction in 
A c possible from IN. 

Figure 1 shows A c accounting for the competition between 
homogeneous and heterogeneous freezing during cloud for- 
mation (Barahona and Nenes, 2009b) (details of the calcu- 
lation are described Sect. 3.3). As A c is strongly influenced 
by the vertical velocity at the point of freezing and much less 
sensitive to small T variations, the curves in Fig. 1 represent 
the peak A c that would be obtained in box models after a sin- 
gle freezing event (Barahona and Nenes, 2008; Hoyle et al., 
2005; Karcher and Lohmann, 2002; Spichtinger and Gierens, 
2009a). If Ain is always very close to An m , competition be- 
tween homogeneous and heterogeneous freezing could yield 
A c close to observations. This requires Ain ~ 0.1 cm -3 , 
which is 20-fold higher than measured dust concentrations 
(~. 005 cm -3 ) at the tropopause level (Froyd et al., 2009). 
Ammonium sulfate aerosol is present at much higher concen- 
trations than dust, and can serve as IN (Abbatt et al., 2006; 
Wise et al., 2010) if a fraction of them is effloresced (which 
is possible, given that it deliquesces at ~90% relative humid- 
ity) (Fountoukis and Nenes, 2007; Shilling et al., 2006). 

To inhibit homogeneous freezing and reproduce obser- 
vations of A c , the concentration of ammonium sulfate IN 
needs to be within 10% of Au m ; if concentrations fall below 
0.9An m , homogeneous freezing is triggered and predicted A c 
is significantly above observations (Fig. 1). If higher concen- 
tration than Aii m is present, homogeneous freezing is com- 
pletely suppressed, but too many crystals still form (Bara- 
hona and Nenes, 2008). In fact, if all ammonium sulfate is 
available as IN, A c from heterogeneous freezing and pure ho- 
mogeneous freezing are always comparable (Fig. 2), because 
crystals formed from ammonium sulfate IN are very small 
(with size close to the dry aerosol; 0.02-0.05 pm, Froyd et 
al., 2009) and grow too slowly to quench supersaturation be- 
fore a large fraction of the aerosol freezes heterogeneously. 
A c is within observed values only if the average size of 
crystals at the point of freezing is 2 pm or larger (Fig. 2), 
which is too large for upper tropospheric aerosol (Froyd et 
al., 2009). Experimental studies suggest that heterogeneous 
freezing of ammonium sulfate IN at T ~ 240 K can be very 
selective (about 1 in 10 5 particles nucleate ice. Shilling et 
al., 2006). If the same selectivity maintains at lower T , too 
few IN would be available to prevent homogeneous freezing 
(therefore resulting in high A c ). Higher nucleation selectivity 
(e.g., about 1 in 10 2 particles actively nucleating ice) would 
result in complete inhibition of homogeneous freezing and 
still maintain A c close to observations (not shown). A pure 



Fig. 1. Ice crystal concentration, A c , as a function of updraft 
velocity, u. The cloud was assumed to form at T = 185 K and 
p = 100 hPa (details provided in Sect. 3.3.1). Low values of u cor- 
respond to cloud formation driven primarily by large scale dynam- 
ics, whereas u > 50 cm s - 1 is characteristic of cirrus developing in 
the vicinity of convective systems with intense gravity wave break- 
ing (Kim et al., 2003). Solid line indicate A c calculated for pure 
homogeneous freezing, dashed line for Ain — A[i m , and dotted for 
Ain — 0.75An m . For Ain — An m , A c lies close to the observed val- 
ues for u < 50cms -1 (Kramer et al., 2009) but is very sensitive to 
small fluctuations in Ain- 

heterogeneous scenario of ice nucleation on ammonium sul- 
fate however implies a maximum supersaturation below 20% 
(Fig. 3b), i.e., S greater than 1.2 would be rarely observed 
as newly formed crystals would rapidly remove supersatura- 
tion. This is at odds with observations of relative humidity 
that suggest that clear-sky supersaturation in the vicinity of 
cirrus up to 70% (and occasionally above) is very common at 
low temperature (Kramer et al., 2009). Hence the low A c and 
high S observed at high-level cirrus can be reconciled with 
box-model results only if the concentration of ammonium 
sulfate IN is remarkable constant (0. 1 ± 0.01 cm -3 ), the con- 
centration of dust is exceptionally large, or, the fluctuations 
in vertical velocity from gravity wave motion are neglected. 

The freezing fraction of organic glassy aerosol can be 
much lower than that of ammonium sulfate and maintain 
Ain close to A i j m (hence yield low A c , Fig. 3a) if the ver- 
tical velocity is below 15cms -1 (Murray et al., 2010). At 
larger updrafts however, homogeneous nucleation is trig- 
gered, producing high A c (Fig. 3a). When integrated over 
a normal distribution of updrafts with standard deviation er u 
(Fig. 3c), A c remains within observed values for a u up to 
40cms -1 at T — 195 K. The onset of homogeneous nucle- 
ation occurs at even lower u for colder temperatures and 
A c deviates from observations for er„ as low as lOcms -1 
at T = 185 K. Predominance of heterogeneous nucleation 
from glassy IN would also imply maximum supersaturation 
around 30% if the new formed ice crystals efficiently remove 
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Fig. 2. Simulations of ice crystal concentration by pure heteroge- 
neous freezing. N c is presented as a function of the initial size of the 
ice nuclei. Conditions (Lawson et al., 2008) used were T = 185 K, 
p = 100 hPa, and vapor-to-ice deposition coefficient, ay, of 0.07 
(dashed line), and 1.0 (solid line). The IN population was assumed 
to be monodisperse with total number concentration of 100 cm -3 
(Lawson et al., 2008). Freezing of solid ammonium sulfate was 
assumed to occur in a “burst” around the heterogeneous freezing 
threshold described by sigmoidal freezing spectrum with inflection 
point Shet = 15% (Abbatt et al., 2006), where 99% of the aerosol 
freeze within a 2% supersaturation interval about Sh et (inset plot). 

supersaturation (Barahona et al., 2010a; Murray et al., 2010). 
However as the freezing fraction of glassy IN is small, S can 
increase even after heterogeneous nucleation has occurred. 
Still, S would have to remain below S^om for N c to remain 
low unless homogeneous freezing is suppressed. The later 
scenario is however not supported by observations. Both, in- 
cloud and clear-sky RH are generally limited by the homoge- 
neous freezing threshold indicating efficient supersaturation 
removal by homogeneously-frozen ice crystals (Kramer et 
al., 2009; Selkirk et al., 2010). Although uncertainty in RH 
can be typically up to 20% (Kramer et al., 2009), it is still 
smaller than the difference between the homogeneous and 
heterogeneous freezing thresholds, typically between 30% 
and 40%, giving this support to the idea that homogeneous 
freezing occurs at low temperatures. All together, this im- 
plies that in the presence of (ubiquitous) T fluctuations, the 
presence of glassy IN may contribute, but not fully account 
for the observed characteristics of low T cirrus. 


3 Parcel statistical ensemble model 



u (cm S ^ 



u (cm s ’) 



ct u (cm s ) 

Fig. 3. Comparison between heterogeneous effects from solid am- 
monium sulfate (Abbatt et al., 2006) and glassy citric acid aerosol 
(Murray et al., 2010), using the analytical model of Barahona and 
Nenes (2009a) for homogeneous and heterogeneous freezing, (a) 
Maximum ice crystal concentration as a function of updraft velocity 
for a single freezing event, (b) Maximum supersaturation achieved 
for a single freezing event, (c) Ice crystal concentration averaged 
over a normal distribution of updraft velocities with zero mean and 
standard deviation <r u . The gray lines represent the range of N c 
typically observed (Kramer et al., 2009). 


The main processes affecting the evolution of N c and mean phase. The magnitude of each process can be expressed in 

saturation ratio, within a cirrus layer are the freezing of terms of a characteristic timescale, i.e., rf r , r se d, tm, and r re i 

new ice, the sedimentation of existing ice crystals, the lift- for freezing, sedimentation, lifting, and relaxation, respec- 

ing of air masses (which generates supersaturation), and the tively. Fluctuations in S and u can have a strong impact on 

relaxation (i.e., mass transfer) of water vapor to/from the ice all cloud processes; we therefore represent them in terms of a 
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probability distribution centered about the cirrus-average sat- 
uration ratio, S 0 , and vertical velocity, Tt. The width of these 
probability distributions is largely determined by the mean 
amplitude of temperature fluctuations, ST (Bacmeister et al., 
1999; Hoyle et al., 2005; Karcher and Burkhardt, 2008). The 
rate of ice production is given by the frequency with which 
S exceeds the homogeneous freezing threshold (Karcher and 
Burkhardt, 2008) times the length and intensity of each freez- 
ing event (hence tf r ) (Barahona and Nenes, 2008; Pruppacher 
and Klett, 1997). The same fluctuations also affect the local 
mass transfer rate between the ice and vapor phases, so that 
when averaged over the cloud, water deposition/sublimation 
occurs at an “effective” saturation ratio, 5 e ff, that may differ 
from S 0 - 

Supersaturation and crystal number in the cirrus cloud 
are determined using a “Lagrangian statistical ensemble” ap- 
proach. This involves determining the time-dependant state 
of i homogeneous Lagrangian “parcels” that move with a 
(time-dependant) vertical velocity, up, ensemble averaging 
of the parcel solutions (outlined below) weighted by the ap- 
propriate probability distribution give approximate equations 
that describe the time-dependant properties for the whole cir- 
rus. From these considerations, simple equations can be de- 
rived that represent the evolution of N c and S 0 in the cirrus 
(Sects. 3.1, 3.2). 

3.1 Evolution of saturation ratio 


In the absence of ice nucleation, the rate of change of satu- 
ration ratio, S, within the 1 th Lagrangian parcel is given by 
(Barahona and Nenes, 2009a; Seinfeld and Pandis, 1998) 


d Si 

= ail; Si 

dt 


OO 

/ 


Dz 


dD c 


dt 


-n c j(D c )dD c 


(1) 


whprr m — S AH s M ™ _ KtZtjL an H v — Pi 71 1,1 '■‘ p A H is thp 

wnere a - ^ RJ , 2 RT and Y - p . d2 m vP °’ ls ttle 

latent heat of sublimation of water, g is the acceleration of 
gravity, c p is the heat capacity of air, p° is the ice satura- 
tion vapor pressure at T (Murphy and Koop, 2005), p is the 
ambient pressure, M w and M a are the molar masses of wa- 
ter and air, respectively, R is the universal gas constant, p, 
and p a are the ice and air densities, respectively, and D c is 
the volume-equivalent diameter of an ice particle (assuming 
spherical shape). n c ,i(D c ) is the ice crystal size distribution 
in the ; th parcel, and 


dD c ,i _ G(Sj- 1 ) 

dt D c ,i 


where G 


PiRT 

4 


A H s pj ( A H S M W 
4k d T ^ RT 



, k a is the 


/ / 

thermal conductivity of air, D v — D v ('f. p.a d \) is the water 
vapor diffusion coefficient from the gas to ice phase corrected 
for non-continuum effects, and ad is the water vapor depo- 


sition coefficient. Substituting Eq. (2) into Eq. (1) provides 
after evaluation of the integral. 


dS, (Si- 1 ) 

—aiij Sj 

dt T re l, 


where r re i j = (yS7V c , D C ,) 1 is the relaxation time scale in 
the i th parcel, f J > — yG, and, N c j, D c j are the concentration 
and mean size of ice crystals in the i tb parcel, respectively. 

Equation (3) provides the supersaturation “state” for every 
Lagrangian parcel considered in the ensemble. Knowledge 
of the distribution of ip (from the spectrum of gravity waves 
in the cirrus) can then be used to “drive” the parcels in the en- 
semble to find the resulting distribution of .S', . Averaging is 
carried out first over all parcels reaching a given cloud level 
with vertical velocity u j (referred to as the “j th cloud veloc- 
ity state”), and then averaging over all cloud states. Based on 
this, the average saturation ratio, S 0 , of the cloud over a time 
interval At is 


So(0 = 


+00 1 

=/ //*<"• 


x , r )P(p,x, r)i irdxdp 


- ooX(l ) 0 


(4) 


where p, = |, u and u are the instantaneous and average verti- 
cal velocity, respectively, x denotes the position in the cloud, 
r = +, where t' is the averaging time, and X(t) is the do- 
main of x. P(p,x, r) is the normalized probability at time 
t' of finding a parcel between position ,v and x+dx (where 
dx — d y dyd " ), with vertical velocity within u and u + du. 

Equation (1) can be simplified, by considering that fluctu- 
ations generated by gravity waves are random in nature (i.e., 
follow a Gaussian distribution. Fig. 4d). Thus, under the as- 
sumption that P(p,x, r) does not vary with space and time 
over At, P(p,x, r) ~ P(p) and Eq. (4) simplifies to 


S 0 (t) = 


+00 1 

/// 


Si(p,x , r)P(p)drdxdp 


-OO X(t) 0 


(5) 


Equation (5) assumes that S Q is affected by processes that act 
throughout the volume of the cirrus cloud. Other processes, 
like entrainment and radiative cooling, are neglected. Al- 
though this will not affect the conclusions of our study, they 
could be included in future studies e.g. indirectly through 
appropriate modification of the vertical velocity distribution 
(Barahona and Nenes, 2007). 

l 

Defining Sj = f J Si(p,x, r)drdx as the average super- 
X(t) o 

saturation of parcels in the “j” velocity state over the time 
interval At, Eq. (5) can be rewritten as 


+oo 

S 0 (t) = J Sj(pj)P(pj)dpj 

— OO 


(6) 


the time derivative of which gives. 
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HI 
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(Si - 1) 

frel,i 


drdx—i 


—auj J J Si 


drdx 


X(t) o 


-/ J — -^jdrdx (10) 

X(t) o 

which can be rewritten as, 

^r = ° Ui ' Sl ~ f 'tH* (11) 

X(t) o 

Introducing ,S' e ff. / so that, 


Fig. 4. Evolution of a cirrus cloud under pure homogeneous freez- 
ing, lifting at 1 cm s~ 1 with initial T = 195 K and cloud thickness, 
H = 500 m, and dj = 1. Shown are (a) the ice crystal number con- 
centration. (b) mean supersaturation, (c) characteristic timescales 
of freezing (gray dots), relaxation (solid lines), and sedimentation 
(dotted lines), and, (d) frequency distribution of vertical velocity, 
for different values of the mean amplitude temperature fluctuations, 
ST. 


dS 0 

dt 


+ 00 _ 
f dSj(fij) 

J 


+oo 

{ 5 , dp 'l*l>. 


— OO 


— CO 


(7) 


the second integral on the right hand side of Eq. (7) de- 
pends on the source of vertical velocity fluctuations. Distant 
sources of gravity waves result in stationary P(/uj), hence 
— * 0. However P(p-j) can be perturbed by near con- 
vective and orographic sources; in such cases P(/zj) is not 
completely Gaussian and exhibits a tail towards high veloc- 
ities (Bacmeister et al., 1999). For the purpose of this study 
it is assumed that dP ^ j) = 0, which implies that the char- 
acteristic amplitude of temperature fluctuations, ST, remains 
constant during the entire period of simulation. Equation (7) 
then becomes 


dS Q 

dt 


+oo _ 



dSj(fij) 

dt 


P(Hj)dnj 


— OO 


( 8 ) 


Using the definition of Sj, 
l 

dSj f f dSi(/j,,x, r) 

dt J J dt 

X(t) o 

Substitution of Eq. (3) into above provides 


X(t) 0 X(t) 0 


1 

drdx 

tf el, / 


( 12 ) 


From Eq. (3), 




X(t) o 


(13) 


Combining Eqs. (12) and (13), Eq. (11) can be written as 


dJl 

dt 


' = otujSj - 


S e ffJ - 1 
r rel J 


(14) 


where r re i j is the relaxation time scale associated with the 
j th state. S e ffj is an “effective” saturation ratio for de- 
position/sublimation processes, defined below. Introducing 
Eq. (14) into Eq. (8), 


+oo , - \ 

-Jf= J l au jSj- ^ 1 j P(jXj)diij (15) 

or, 

+oo +oo / _ V 

J otujSjPinjWiXj- J ( ^ . 1 ) p( Pi) d Pj ( 16 ) 

-OO -OO ' Te ’ J / 


The first term in the right hand side of Eq. (16) must be equal 
to u S 0 , as in the absence of deposition/sublimation, S 0 in 
the layer increases exponentially with time (Pruppacher and 
Klett, 1997). With this, Eq. (16) becomes. 


d S 0 
dt 


So 


Bift 



SeffJ ~ 1 

tfelj 


P(pj)dtLj 


(17) 


where t| , f t = (au)~ 1 . Equation (17) must be solved for each 
time step specifying P (fi ,■) and then evaluating S e ff, ; and the 
integral on the right hand side. Since P(/J,j) is determined 
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by the random overlapping of gravity waves of different fre- 
quency and amplitude (e.g., uj is given by a Fourier series 
in time. Sect. 3.3), then for a time step of integration much 
smaller than ri;f t (~10 2 s) Eq. (17) can be approximated by 


dS 0 5eff j 

—— — an 8 o — 

at T rc \j 


(18) 


U=Uj 


where 5 e ff, / and r re i, / are calculated at the instantaneous ver- 
tical velocity. 

Equation (18) gives the evolution of the S 0 in the cirrus 
cloud; its solution however requires the knowledge of 5 e ff, / • 
This is accomplished by considering the properties of the dif- 
ferent parcels reaching the cloud layer at t. For example, if 
Si in the i th parcel is a pseudo-steady state, ^ ~ 0 (Korolev 
and Mazin, 2003) and from Eq. (3), 


Si , ss — 


Bift ,i 


flift.i Ael 


(19) 


where run,; = (a w,) -1 and S), ss is the steady state satura- 
tion ratio in the I th parcel. If rm, , < 0 then .S', ss < 1, and 
vice-versa. Thus, if u j < 0, the layer would likely be sub- 
saturated over At (e.g., Eq. 6), and vice-versa when u ; > 
0. Thus, depending on the sign of uj there is net deposi- 
tion/sublimation of water vapor in the cloud layer. Not all 
parcels however reach steady state; therefore the degree of 
saturation/subsaturation associated with the j tb state depends 
on the probability distribution of saturation within the cloudy 
layer, P s (S,S 0 ,ST), which is a function of S 0 and the aver- 
age amplitude of temperature fluctuations, ST. Thus, ,S’ C n for 
Uj <0 is found by averaging over all states that would lead 
to subsaturation, i.e., P s (S,ST,S 0 ) for which S < 1. Sim- 
ilarly, when uj >0, the supersaturated (S > 1) region of 
P s (S,8T , S 0 ) is used. 


fs 


dP s jS,ST,S 0 ) 
dS 


dS 


Sc IT. j — 


( 20 ) 


j dP s <S,ST. So ) ds 


where a= 


1 u j > 0 

0 Uj < 0 


and b — 


3'hom tt j > 0 
1 Uj < 0 


The homogeneous freezing threshold, .S'| lom , is set as the up- 
per limit of P s (<57\ S 0 ) as ice crystal production quickly re- 
moves supersaturation above Shorn (Karcher and Burkhardt, 
2008; Karcher and Haag, 2004). 

3.2 Evolution of ice crystal number concentration 

The evolution of the number concentration within a cloudy 
layer is given by 


dN c _ dN c 
dt dt 


dNc 

dt 


( 21 ) 


where 


the layer, and 


is the rate production of ice crystals within 


sed 


is their sedimentation rate. If homo- 


geneous and heterogeneous nucleation are active, 
given by. 


dNc 

dt 


_ dNc 

fr dt 


dN c 

fr.hom 


where ^ 


fr,hom 




fr,het 

are the ice crystal produc- 


fr,het 


tion rates from homogeneous and heterogeneous nucleation, 
respectively. 

3.2.1 Homogeneous nucleation 

Ice crystal production by homogeneous nucleation is driven 
by local motions and occurs within single parcels when 
Si > Shorn- The maximum ice crystal concentration produced 
by homogeneous freezing within the i th parcel is given by 
(Barahona and Nenes, 2008; Pruppacher and Klett, 1997) 


Ncj = N 0 


1 — exp 


( ^max. 

-/ 

o 


v 0 J (Sj)dt 


( 22 ) 


where f max ,; is the time at which ice crystal nucleation stops, 
J is the homogeneous nucleation rate coefficient and N 0 , v 0 
are the deliquesced aerosol number concentration and av- 
erage volume, respectively. Taking the time derivative of 
Eq. (22) gives. 


dN, 


C ,1 


dt 


= N 0 v 0 J (Si)exp 


fr,hom 



v 0 J(Sj)dt 


(23) 


which can be approximated by (Barahona and Nenes 2008) 


dN, 


C ,1 


dt 


' NoVoJ max, i ex P 


fr,hom 



max.i = 7(S m ax,!'). S max ,! is the maximum saturation 


where ./, 

ratio reached in the ; th parcel, calculated by setting — 0 
inEq.(l), 


Siri ax./ — 


Y_ 

au 


CX. 

-/ 


Dt 


dD c 


dt 


n i, nuc (Dc)dDc 


(25) 


sed 


where «/, nuc (D c ) is the size distribution of the recently nucle- 
ated ice crystals, and D 0 is the mean size of the deliquesced 
aerosol. Equation (25) assumes that only recently nucleated 
ice crystals are contained within the parcel. In reality, a frac- 
tion of preexisting crystals remain in nucleation zones (typ- 
ically located near the cloud top, Spichtinger and Gierens, 
2009b) inhibiting the homogeneous freezing of ice. Ice crys- 
tals experience gravitational settling, hence only those crys- 
tals with terminal velocity, « te rm. below u would be found at 


www.atmos-chem-phys. net/1 1/3757/201 1/ 


Atmos. Chem. Phys., 11, 3757-3771, 2011 



3764 


D. Barahona and A. Nenes: Dynamical states of low temperature cirrus 


the cloud top. Adding the consumption of water vapor from 
preexisting crystals to the right hand side of Eq. (25) gives 


S ■ — ^ 

aui 


CX. 

/ 

LP 0 


q d Dr i 

D- ci —^ni, auc (D c )dD c 


Dterm 

[ 9 dD c 

+ / D 2 c -^n c (D c )dD c 


Drr 


(26) 


where n c (D c ) is the cloud ice crystal size distribution, D teim 
is the size of the crystal for which M term = u, and D m ; n is 
the minimum size of the preexisting crystals in the cloud. 
Equation (26) can be combined with Eq. (2) to obtain 


S — ^ 

Jmax.i — 


LXJ 

i Dl 


dD c 

dt 


d nuc ( D c . .S'niax./ )d D c 


+ GA c D c / ps (5 m ax,i-l)] 


(27) 


where, / ps : 


Dterm 


N C D C 


D c n c (D c )dD c , is the fraction of pre- 


existing ice crystals remaining in nucleation zones. As ice 
crystals remaining in the cloud layer were produced by pre- 
existing freezing events, Eq. (27) provides a link between the 
history of different parcels and the nucleation of new crystals. 
The analytical solution of Eq. (27) is presented elsewhere 
(Barahona and Nenes, 2009b; Barahona et al., 2010b). 

The rate of ice crystal production by homogeneous nucle- 
ation in the j th cloud velocity state is given by the concentra- 
tion of nucleated crystals over the freezing timescale. 


dN, 




dt 


where 


fr.hom 


Ps(S > Shorn) N 0 

^hom, j 


Hy (U j ) 


(28) 


S’ max 


hom. 


j — U 0 /max,; exp I f J(Sj)dSj 


H v ( u j ) is the Heaviside function and is introduced because 
homogeneous nucleation is very unlikely in parcels with 
negative vertical velocity (i.e., updraft must be maintained 
for some time before Shorn is reached after which it is quickly 
depleted by crystal nucleation and growth, Barahona and 
Nenes, 2008; Karcher and Lohmann, 2002) . P s (S > Shorn) 
represents the fraction of parcels for which S > Shorn- Using 
the same averaging procedure as for the supersaturation 
equation, we obtain 


dN c 

dt 


= No Ps(S> Sh om ) ^ 

fr.hom r homJ 


(29) 




3.2.2 Heterogeneous nucleation 


The formulation of the ice crystal production rate by het- 
erogeneous freezing is simplified by using the ice nucle- 
ation spectrum, /Vin(S, T) (Barahona and Nenes, 2009a) . 


Assuming that N\n(S,T) is weakly dependent on T and 
/Vin(S, T ) = fV 0 ,het/het(S), then following Eq. (29) we write. 


dN c 

dt 


fr,het 


d /het d S i 

= AUet^ — / J s(S > S h et)tfv(Mj)|^ (30) 


dS dt 


where, iV 0 ,het is the total number concentration of IN species, 
and /het, Shet its freezing fraction and heterogeneous nucle- 
ation threshold, respectively. As Ain is usually small, S is 
not immediately depleted by ice crystal growth and ^ can 
be approximated by the instantaneous rate of increase of su- 
persaturation (this is further justified in the case of glassy IN 


as 


dfba (S) ; 


dS 


dN c 

dt 


is constant , Murray et al., 2010). Therefore, 

=N oMx d ^ocujS P s (S>S het )H v (^j)\ (31) 
dS 'd—d-j 


Ifr.het 

which can be written as 
dN 


dt 


Hy (LL j ) 

= iVo,hetA(S>Shet)^^ 

fr.het r het,j 

, -1 


(32) 




where thet J = (^jir au jS) ■ 

3.2.3 Sedimentation of ice crystals 

Sedimentation processes out of the cloud layer depend pri- 
marily on the bulk properties of the cloud, i.e., the mean ice 
crystal size distribution and number concentration (interac- 
tion of individual parcels with falling crystals within the layer 
is accounted for in Eq. 27). The ice crystal loss rate by sedi- 
mentation is then given by. 


dN c 

dt 


sed 


oc 

-hf 

£>min 


— TT I Mterm(77c)n(D c )f/D c 


(33) 


where H is the cloud layer thickness. As Mterm ~ Dc (Heyms- 
field and Iaquinta, 2000), Eq. (33) can be further simplified 
to 


dN c 

dt 


sed 


-^c^term 

H 


Nc_ 

^sed 


(34) 


where Mterm — M term (77c)- 

3.3 Numerical solution 

3.3.1 Competition between homogeneous and 
heterogeneous freezing 

Calculation of ice crystal number concentration, N c in in- 
situ cirrus from combined homogeneous and heterogeneous 
nucleation in Figs. 1 and 3 is done using an analytical pa- 
rameterization developed for in situ formed cirrus clouds and 
freezing fractions below 0.6 (Barahona and Nenes, 2009b). 
When the calculated freezing fraction exceeds 0.6, a sig- 
moidal increase in N c is assumed (Barahona et al., 2010a), 
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in agreement with parcel model simulations and field obser- 
vations (Barahona and Nenes, 2008; DeMott et al., 2003). 
For combined homogeneous and heterogeneous nucleation, 
it was assumed that the IN freeze instantaneously at a super- 
saturation freezing threshold, .Vh c t, of 15%, typical of deposi- 
tion mode IN (Abbatt et al., 2006), with a 0.1 pm diameter at 
freezing (Froyd et al., 2009). Glassy aerosol was assumed to 
have a total concentration of 50 cm -3 and a freezing fraction 
given by the nucleation spectrum of Murray et al. (2010). 

3.3.2 Vertical velocity spectrum 

Observations suggest that the spectrum of gravity wave per- 
turbations near the tropopause is pseudorandom in nature, a 
result of the superposition of waves from different sources, 
and therefore varies temporally and spatially (Bacmeister et 
al., 1999; Jensen and Pfister, 2004; Kim et al., 2003; Sato, 
1990). On average, the associated vertical velocity spec- 
trum can be approximated by a Gaussian function, although 
this may underestimate the frequency of high amplitude per- 
turbations (Bacmeister et al., 1999; Herzog and Vial, 2001; 
Kim et al., 2003). A representative spectrum of vertical ve- 
locity fluctuations can be generated using a Fourier series 
(Bacmeister et al., 1999; Jensen and Pfister, 2004) of the 
form u = u + '%2A(ztTj)cos(Tnjt + mH + cp) where m is the 
j 

vertical wave number, H is the cloud thickness, and tuj, 
A(uij), and cp, are the wave frequency, phase, and amplitude, 
respectively. We have adopted this representation as follows. 
For each cirrus simulation, a time series of u was generated 
over the frequency interval uj = [3.35 x 10 7 .9.44 x 1 0 4 ] 
Hz (Jensen and Pfister, 2004), using randomly generated <p 
and m. A(nr ; ) was calculated using a power spectrum scal- 
ing law of —1.85 for uij > lx 10~ 5 Hz and of —0.25 for 
m j < 1 x 10~ 5 Hz (Herzog and Vial, 2001; Jensen and Pfis- 
ter, 2004). This procedure resulted in a normal distribution 
of u (Fig. 4d) centered around u. The maximum amplitude 
was assumed to occur at nij = lx 10 -3 Hz (Jensen and Pfis- 
ter, 2004) as it reproduces the results of Gayet et al. (2004) 
(Fig. 4 green line) which give positive u around 0.23 ms -1 
for ST — 1 K (i.e., A(1 x 10 -3 ) 2AST). Representative 
time series for u(t ) are presented in Fig. 5a. 

3.3.3 Ice crystal production 

The homogeneous freezing timescale, thom. /, was calculated 
using the parameterization of Barahona and Nenes (2008, 
2009b, a). Precursor aerosol was assumed to be composed of 
deliquesced ammonium sulfate, lognormally distributed with 
dry mean geometric diameter of 40 nm, geometric dispersion 
of 2.3, and number concentration of 100 cm~ 3 (Lawson et 
al., 2008). To account for possible compositional impacts 
on crystal growth kinetics, the water-vapor deposition coef- 
ficient was varied between 0.006 (Magee et al., 2006) and 
1.0. Homogeneous freezing is described using the param- 
eterization of Koop et al. (2000). The term P s (S > Shorn) 



Fig. 5. Time series of (a) updraft velocity, u, (b) total water con- 
tent, qt 0 t, and D c (c) mean ice crystal diameter, D c , for the condi- 
tions presented in Fig. 4; (d) time series of D c for initial S and N c 
of —0.2 and 0, respectively, and homogeneous and heterogeneous 
nucleation active (conditions similar to Fig. 8c). 


in Eq. (28) is the probability of finding S above Shorn, and 
is introduced to account for the threshold behavior of ho- 
mogeneous freezing (Karcher and Burkhardt, 2008; Koop et 
al., 2000). The effect of preexisting ice crystals on freez- 
ing was accounted for by allowing a fraction of N c to de- 
plete water vapor and increase rhomj (Barahona and Nenes, 
2009b; Barahona et al., 2010b). The fraction of preexist- 
ing crystals remaining in freezing zones was calculated as 

Altim 

,/ps = -TTjy f D c n c (D c )d D c , where n(D c ) is the ice crys- 

C C Anm 

tal size distribution, D rnnl is the minimum pre-existing crys- 
tal size, and Aerm is the crystal size for which its termi- 
nal velocity, « t erm, is equal to the uplift velocity of the 
cirrus later, ii. « t erm was calculated assuming ice crystals 
have columnar shape with maximum dimension equal to D c 
(Heymsfield and Iaquinta, 2000). Following Heymsfield and 
Platt (1984) it was assumed n(D c ) = AD" 315 ; the parame- 
ters A and D m ; n were calculated from the moments of n(D c ): 
N c = /“ n( D c )d D c and D c = 4~ fn D c n(D c )dD c . The 
calculation of D c is described below. Integration of equations 
of Eqs. (18) and (21) was accomplished using a fixed time 
step of 2 s. Initial values for N c , o = 0.01 cm -3 and So = 1.0 
were set. Sensitivity to using different initial values affected 
only the time required to establish dynamic equilibrium (by 
a few hours) and is assessed in Fig. 8. 

The heterogeneous freezing timescale, Thetj, (e.g., Eq. 32) 
was calculated at each time step using the instantaneous ver- 
tical velocity, uj , and the effective supersaturation S e ff, /. For 
the simulations presented in Sect. 4, the aerosol freezing frac- 
tion, /het, was calculated from the spectrum of Murray et 
al. (2010). N 0 het and 5h et were set to 1 cm -3 and 1.35, re- 
spectively (Murray et al., 2010). 
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3.3.4 Ice crystal sedimentation 


The rate of ice crystal sedimentation over the cloud scale, 
//, was assumed proportional to the terminal velocity of the 
mean crystal size D c (Eq. 33). Other removal processes (ice 
crystal sublimation and detrainment) are neglected; H how- 
ever was varied over a wide interval (100 to 5000 m) to ac- 
count for the uncertainty associated with neglecting these 
processes. D c was calculated so that the total water vapor 
in the layer was partitioned between ice and vapor phases, 

i n °S M 

where q ice = q tot - ' " w , p, is the 


/ 6gjge \ 
\npiN c ) 


i.e., D c — 

ice density (Pruppacher and Klett, 1997), R is the univer- 
sal gas constant, M w is the molecular mass of water, and 
p° is the saturation water vapor pressure over ice (Murphy 
and Koop, 2005); the minimum ice crystal size was set to 
4 pm in agreement with theoretical studies and experimen- 
tal observations (Barahona and Nenes, 2008; Durran et al., 
2009; Kramer et al., 2009). Loss of total water content, 
</tot. from the cloudy layer is also accounted for by solution 
of — — \pi D\ ■ Representative time profiles of 

sed 


Dc and // tot are presented in Fig. 5. The timescale of relax- 
ation at u — iij, t re i. j, was calculated using N c and D c of the 
cloud layer (e.q., Eq. 13). Mass transfer limitation from non- 
continuum effects are taken into account in the calculation of 
Thom, / but neglected in the calculation of r re i j. The latter is 
justified as D c > 10 pm. 


4 Cirrus in dynamical equilibrium 

The model developed in Sect. 3 is used to describe the evolu- 
tion of a cirrus layer at low temperature taking into account 
the effect of internal S variations on N c . Forward integration 
of Eqs. (18) and (29) and (34) is carried out using the proce- 
dure described in Sect. 3.3. A wide range of initial conditions 
and model parameters are selected to describe the cirrus evo- 
lution under different scenarios. 

Figure 4 shows the evolution of a cirrus layer subject to 
gravity-wave fluctuations with an initial average temperature 
of 195 K and lifting at u — lcms -1 . Only homogeneous 
freezing is considered (heterogeneous freezing is “switched 
off”, i.e., fVo het = 0). For values of <5 T > 1 K, the cloud ini- 
tially experiences a strong homogeneous nucleation pulse, 
so that N c initially increases steeply (Fig. 4a); the consump- 
tion of water vapor by crystal growth decreases S 0 (Fig. 4b) 
which prevents any new freezing events. N c slowly decreases 
from sedimentation loss; only after enough ice crystals sed- 
iment out of the cloud layer, S 0 increases (e.g., ST — 1 K, 
green lines) and new freezing events occur. For ST > 1 .4 K 
(purple lines) this is possible even if the layer remains on 
average slightly subsaturated (So ~ 1) because the probabil- 
ity distribution of S is broad enough for a non-negligible 
probability with S > Shorn- However it is likely that recently 
formed crystals will sublimate within a few hours in the sub- 
saturated environment returning the moisture to the layer (ice 


crystal removal by sublimation is not considered). The cir- 
rus is then maintained by new, independent freezing events. 
This “pulse-decay” behavior is characterized by r se d » r re i 
so ice crystals reside long enough in the cloud to relax super- 
saturation (Fig. 4c); this behavior is also consistent with the 
parcel model concept of cirrus (where high N c and low S 0 
coexist within the parcel). The subsaturation levels (Fig. 4b) 
achieved in this state are in agreement with in situ observa- 
tions of relative humidity in dissipating clouds (Gao et al., 
2004; Kramer et al., 2009). 

The cirrus evolution is however quite different when 52” is 
small; the distribution of S is narrow, and substantial ice pro- 
duction is only possible after enough supersaturation (i.e., 
S 0 ) builds up in the cloudy layer to allow a non-negligible 
probability where S > Shorn- Thus, freezing events produc- 
ing large N c (associated with large u fluctuations; Fig. 4d) are 
less frequent. Low N c allows the formation of large ice crys- 
tals (Fig. 5c, d) which sediment out of the layer before sub- 
stantially depleting supersaturation, leading to new freezing 
events. This “dynamic equilibrium” between ice production 
and loss is a previously unidentified microphysical regime of 
cirrus, characterized by r se d ~ r re i j (Fig. 4c); it maintains 
low N c and high S Q in the cloudy layer (Fig. 4a, b) and 
is consistent with observations of low-temperature cirrus. 
Clouds in “dynamic equilibrium” also exhibit broad crys- 
tal size distribution, because large ice crystals coexist with 
freshly-formed (small) crystals in the cloud. The average 
ice crystal size in this case converges to values around 15 to 
20 pm, in agreement with observations (Kramer et al., 2009). 
Before equilibrium is reached, D c exhibits larger values due 
to the model’s sensitivity to initial conditions (Fig. 5c, d). 

When simulations (such as those of Fig. 4) are placed on a 
“state diagram” of N c vs. S 0 , the two microphysical regimes 
described above clearly emerge. Examples are presented in 
Figs. 6 and 7 for a range model parameters and a variety 
of 5 2” (lines of distinct color). For example, decreasing the 
cloud thickness to 100 m (Fig. 4b) increases the sedimenta- 
tion rate (Eq. 34) allowing S to replenish quickly and facil- 
itating the onset of equilibrium. Similarly, the rate of ice 
crystal growth increases with T (Fig. 4d) increasing the av- 
erage ice crystal size (hence decreasing r re i and r se d) and fa- 
cilitating the onset of equilibrium. Dynamical equilibrium 
would however lead to lower N c than observed at 7’ ~ 225 K 
(Gayet et al., 2004; Kramer et al., 2009) and it is likely that 
clouds at these conditions would evolve following a pulse- 
decay behavior. Figure 7 shows state diagrams for different 
values of w, T and H for the similar initial conditions. In gen- 
eral, progression towards a “dynamic equilibrium” is favored 
when supersaturation replenishes quickly (i.e., at high it) as 
ice crystal growth and sedimentation are favored (leading to 
low r se d), and vice-versa for “pulse-decay” behavior. Rapid 
consumption of supersaturation by the growing ice crystals 
also decreases the time between freezing pulses and replen- 
ishment of S. Thus the period of the S and N c oscillations in 
the dynamic equilibrium state is mainly controlled by r re i. 
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Fig. 6. Sensitivity of N c and S 0 evolution to cloud formation con- 
ditions for different values of ST (color scheme same as in Fig. 4); 
(a) same conditions as in Fig. 4, (b) cloud thickness, H = 100 m 
(increased ice crystal removal rate), (c) deposition coefficient equal 
to 0.006 (Magee et al., 2006) (slow water vapor transfer), and (d) 
initial temperature 225 K and cloud lifting at 5 cm s _1 . The yellow 
star in each panel indicates initial conditions. The arrows indicate 
the temporal progression along each trajectory. The integration time 
was 40 h cases, except in (d) were it was 15 h. 


Clouds in the “dynamic equilibrium” regime are also much 
less sensitive to slow water vapor deposition than predicted 
by box-model simulations. Figure 6c shows that dynamic 
equilibrium is possible even for « ( | as low as 0.006 (com- 
pared to ctu = 1 used in Figs. 4 and 6a for the same sim- 
ulation conditions). Still, the high rate of production of ice 
crystals for a t i = 0.006 increases N c and decreases D c , there- 
fore increasing r se d and slowing the replenishment of super- 
saturation. Thus, dynamic equilibrium is only possible for 
ST < 0.8 as opposed to ST <1.0 for aa = 1. Figure 6c how- 
ever shows that the existence of strong kinetic limitations to 
the diffusional growth of ice crystals cannot be ruled out. 

4.1 Effect of heterogeneous IN and initial conditions 

It is important to study the sensitivity of the dynamical states 
of cirrus to initial conditions used in the simulation and to the 
presence of heterogeneous IN. Figure 8a shows S and N c for 
the conditions of Figs. 4 and 6a but starting at S 0 j n i — —0.4 
and /V c o = 0. For this set of initial conditions, the onset of 
oscillating behavior is delayed by a few hours before super- 
saturation is reached. This implies that the temperature of the 
first freezing pulse is lower than in the case with ,S ojrn = 1, 
slightly increasing its strength; the system however dampens 
out these variations and eventually follows a similar trajec- 
tory as in Fig. 6a. 


’] r 
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Fig. 7. Similar to Fig. 6, but varying cloud mean vertical velocity, u, 
initial layer temperature, T a , cloud thickness, H, and mean ice crys- 
tal terminal velocity, t/term- The yellow star in each plot indicates 
initial conditions. The integration time was 40 h in for u= 1 cm s _ 1 
and 1 5 h for u. = 5 cm s~ 1 . 


Figure 8b and c present simulations where both homoge- 
neous and heterogeneous nucleation are active for different 
initial conditions. IN are assumed to originate from glassy 
aerosol with N 0i h e t = 1 cm -3 . By allowing ice crystal for- 
mation at low S, heterogeneous IN decrease the ice crystal 
production rate by homogeneous nucleation hence lower the 
maximum N c . This makes the system more stable to high 
amplitude vertical velocity perturbations, facilitating the on- 
set of dynamic equilibrium, which can be maintained up to 
ST ~ 1.2. The “stabilizing” effect of IN is exemplified in 
Fig. 8c where homogeneous freezing was suppressed, i.e., 
only heterogeneous freezing is active. This leads to a damp- 
ened response to vertical velocity fluctuations so that dynam- 
ical equilibrium is possible even for ST as high as 1.8 K. 
Thus, in a cirrus cloud where heterogeneous nucleation is 
dominant, dynamical equilibrium is very robust against ex- 
ternal perturbations. The values of N c and S at the equi- 
librium state are however not significantly influenced by the 
initial conditions nor the ice nucleation mechanism, e.g., af- 
ter sometime the system oscillates about the same values as 
in Figs. 6 and 7. This means that the presence of IN can 
“help” the cirrus to achieve dynamical equilibrium without 
modifying the equilibrium values. 
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Fig. 8. Similar to Fig. 6, varying initial S and N c and the active 
freezing mechanism. 

Figures 6 to 8 also show that the “dynamic equilibrium” 
state occurs spontaneously when ST goes below a charac- 
teristic transition value (which depends on u, T and the pre- 
dominant freezing mechanism). It can also be reached after a 
cloud initially resides in a “pulse-decay” state, if <5T is close 
to the characteristic value (ST ~ 1 K in Fig. 6). When max- 
imum N c and time-averaged S 0 are presented on the state 
diagram for all simulations considered where homogeneous 
freezing is active, the conditions of ST that separate “pulse- 
decay” and “dynamic equilibrium” regimes seem to be uni- 
versal (Fig. 9). If only heterogeneous nucleation is active, 
the cloud resides mostly in the dynamic equilibrium state. 

5 Conclusions and implications 

From the discussion above, cold cirrus clouds will reside in 
the “dynamic equilibrium” regime if ST is below a charac- 
teristic threshold. High-amplitude, orographically-generated 
gravity waves are ubiquitous (Kim et al., 2003) but often 
lose intensity with altitude, weakening their contribution to 
the background spectrum of temperature fluctuations. ST 
can thus decrease enough at high altitude for cirrus to transi- 
tion from a “pulse-decay” to a “dynamic equilibrium” state 
(Fig. 9). This would explain why low N c and high S 0 are ob- 
served at low temperatures near the tropopause. Dynamical 
equilibrium is also possible at warmer conditions (particu- 
larly for high u; Fig. 6d) but require small ST; given that 
high amplitude fluctuations are widespread at lower altitudes 
(Hoyle et al., 2005), cirrus clouds are likely forced to always 
follow a pulse-decaying behavior. Heterogeneous IN never- 
theless may help to stabilize the system so that the dynamic 
equilibrium manifests at higher ST than for clouds with pure 
homogeneous freezing. 



N c (cm 3 ) 

Fig. 9. Maximum ice crystal concentration obtained during the 
cloud evolution simulations against the time-averaged mean satu- 
ration ratio. Results presented for all simulations carried out in this 
study. Integration time varied between 15 and 40 h. Symbols are 
colored by the value of ST used. Regions where the cloud sponta- 
neously transitions to a “pulse-decay” and “dynamic equilibrium” 
state are noted; the “transitional” region marks where the cloud gen- 
erally initially exhibited “pulse-decay” behavior over few hours and 
then transitioned to a “dynamic equilibrium” regime. (O) Homo- 
geneous freezing is active; (A) Homogeneous freezing suppressed. 


In summary, cirrus clouds at low temperature exhibit char- 
acteristics (e.g., low N c and sustained high saturation ratios) 
that cannot be explained with the simple “conventional” pic- 
ture of homogeneous freezing driven by expansion cooling 
in the presence of ubiquitous temperature fluctuations. Even 
if heterogeneous nucleation is dominant, conventional mod- 
els of cirrus could explain the characteristics of low-T cir- 
rus only for weak updrafts, and require neglecting the higher 
amplitude components of the vertical velocity spectrum. In- 
stead, we show that small-scale fluctuations from the action 
of gravity waves can switch a cloud into a previously un- 
known “dynamic equilibrium” regime, with sustained levels 
of low N c and high saturation ratios consistent with “puz- 
zling” characteristics observed in low temperature cirrus. 

With this study, a new understanding for cirrus clouds 
emerges, where the “unperturbed” microphysical state is one 
of dynamical equilibrium with low crystal number and high 
supersaturation. Only when the mean amplitude of tem- 
perature fluctuations exceeds a threshold value (ST ~ IK 
when homogeneous freezing is active) cirrus exhibit the well- 
known “pulse-decay” microphysical state. Throughout much 
of the atmosphere, the latter state dominates, simply because 
ST is larger than the characteristic threshold value. In the 
TTL, ST is still remarkably large (0.6-0. 8 K) (Bacmeister et 
al., 1999; Jensen and Pfister, 2004; Sato, 1990), but does 
not systematically exceed the threshold for “pulse-decay” 
behavior, so cirrus regress to their “unperturbed” dynamic- 
equilibrium state. The presence of heterogeneous IN can 
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also dampen the effect of vertical velocity perturbations on 
ice crystal production facilitating the transition to dynamic 
equilibrium. However, the existence of IN is not a neces- 
sary condition to explain the characteristics of cirrus clouds 
at low temperature. Figure 6c shows that dynamic equilib- 
rium states can also exist even under conditions of strong ki- 
netic limitations (very low ad) and therefore their existence 
cannot be ruled out based on simplified models of cirrus for- 
mation. 

The structure and responses of cirrus to dynamical and mi- 
crophysical forcings can also be portrayed. For example, 
cirrus formed in the region of convective anvils might ex- 
hibit “pulse-decay” state until gravity-wave fluctuations de- 
cay to below the ST threshold and transition to a dynamic- 
equilibrium state. For the same reasons, IN impacts on cirrus 
properties can be strong for clouds in pulse-decay state, but 
not for clouds in dynamic equilibrium; e.g., IN can force the 
cloud to fall towards equilibrium however will not modify 
the equilibrium state. In conclusion, the discovery of dy- 
namic equilibrium states reshapes our understanding of cir- 
rus clouds and their role in anthropogenic climate change, as 
the type of dynamical forcing and the presence of IN will set 
these clouds in one of two “preferred” microphysical regimes 
with very different susceptibility to anthropogenic aerosol. 

Acknowledgements. This study was supported by NASA ACMAP 
and a NSF CAREER award. 

Edited by: J. H. Seinfeld 

References 

Abbatt, J. P. D., Benz, S., Cziczo, D. J., Kanji, Z., and Mohler, O.: 
Solid ammonium sulfate as ice nuclei: a pathway for cirrus cloud 
formation, Science, 313, 1770-1773, 2006. 

Bacmeister, J., Eckermann, S. D., Tsias, A., Carslaw, K. S., and 
Peter, T.: Mesoscale temperature fluctuations induced by a spec- 
trum of gravity waves: A comparison of parameterizations and 
their impact on stratospheric microphysics, J. Atmos. Sci., 56, 
1913-1924, 1999. 

Barahona, D. and Nenes, A.: Parameterization of cloud 

droplet formation in large scale models: including ef- 

fects of entrainment, J. Geophys. Res., 112, D16026; 
doi: 16010. 1 1029/16207JD008473, 2007. 

Barahona, D. and Nenes, A.: Parameterization of cirrus formation 
in large scale models: Homogeneous nucleation, J. Geophys. 
Res., 1 1 3, D 1 1 2 1 1 , doi: 1 1 2 10. 1 1029/1 2007JD009355, 2008. 
Barahona, D. and Nenes, A.: Parameterizing the competition be- 
tween homogeneous and heterogeneous freezing in ice cloud for- 
mation polydisperse ice nuclei, Atmos. Chem. Phys., 9, 5933- 
5948, doi: 10.5 194/acp-9-5933-2009, 2009a. 

Barahona, D. and Nenes, A.: Parameterizing the competition be- 
tween homogeneous and heterogeneous freezing in cirrus cloud 
formation - monodisperse ice nuclei, Atmos. Chem. Phys., 9, 
369-381, doi: 10.5 194/acp-9-369-2009, 2009b. 

Barahona, D„ Rodriguez. J. M., and Nenes, A.: Sensitivity 

of the global distribution of cirrus Ice crystal concentration 


to heterogeneous freezing, J. Geophys. Res., 115, D23213, 
doi:232 10. 2 1029/220 10JD0 14273, 2010a. 

Barahona, D., West, R. E. L., Stier. P„ Romakkaniemi, S., Kokkola, 
H., and Nenes, A.: Comprehensively accounting for the effect of 
giant CCN in cloud activation parameterizations, Atmos. Chem. 
Phys., 10, 2467-2473, doi: 10.5 194/acp- 10-2467-2010, 2010b. 

DeMott, P. J., Meyers, M. P, and Cotton, R. W.: Parameteriza- 
tion and impact of ice initiation processes relevant to numerical 
model simulations of cirrus clouds, J. Atmos. Sci., 51, 77-90, 
1994. 

DeMott, P. J., Cziczo, D. J., Prenni, A. J., Murphy, D. M., Krei- 
denweis, S. M., Thompson, D. S., Borys, R., and Rogers, D. C.: 
Measurements of the concentration and composition of nuclei for 
cirrus formation, Proc. Natl. Acad. Sci. USA, 100, 14655-14660, 
2003. 

Durian, D. R.. Dinh, T.. Ammerman, M., and Ackerman, T.: The 
mesoscale dynamics of thin tropical tropopause cirrus. J. Atmos. 
Sci., 66, 2859-2873, doi:10.1175/2009JAS3046.1, 2009. 

Fountoukis, C. and Nenes, A.: ISORROPIA II: a com- 

putationally efficient thermodynamic equilibrium model for 
K+Ca 2 + Mg 2 +NH+Na + S 04 _ N 0 “CrH 2 0 aerosols, Atmos. 
Chem. Phys., 7, 4639^1659, doi: 10.5 194/acp-7-4639-2007, 
2007. 

Froyd, K. D., Murphy, D. M., Sanford, T. J., Thomson, D. S., Wil- 
son, J. C., Pfister, L., and Lait, L.: Aerosol composition of the 
tropical upper troposphere. Atmos. Chem. Phys., 9, 4363 — 4385, 
doi: 10.5 194/acp-9-4363-2009, 2009. 

Gao, R. S., Popp, P. J., Fahey, D. W., Marcy, T. P, Herman, R. L., 
Weinstock, E. M.. Baumgardner, D. G., Garrett, T. J.. Rosenlof, 
K. H„ Thompson, T. L„ Bui, P. T„ Ridley, B. A., Wofsy, S. C., 
Toon, O. B., Tolbert, M. A., Karcher, B., Peter, T., Hudson, P. 
K., Weinheimer, A. J., and Heymsfield, A. J.: Evidence that nitric 
acid increases relative humidity in low-temperature cirrus clouds. 
Science, 303, 516-520, 2004. 

Gayet, J. F., Ovarlez, J., Shcherbakov, V., Strom, J., Schumann, 
U., Minikin, A., Auriol, F., Petzold, A., and Monier, M.: Cirrus 
cloud microphysical and optical properties at southern and north- 
ern midlatitudes during the INCA experiment. J. Geophys. Res., 
109, D20206, doi: 202 10.21 029/22004 JD004803, 2004. 

Gensch, I. V., Bunz, H., Baumgardner, D. G., Christensen, L. E.. 
Fahey, D. W„ Herman, R. L., Popp, P. J., Smith, J. B., Troy, R. 
F., Webster, C. R., Weinstock, E. M., Wilson, J., Peter, T., and 
Kramer, M.: Supersaturations, microphysics and nitric acid par- 
titioning in a cold cirrus cloud observed during CR-AVE 2006: 
an observation-modelling intercomparison study, Environ. Res. 
Lett., 3, 035003, doi: 10. 1088/1748-9326/3/3/035003, 2008. 

Haag, W., Karcher, B., Strom, J., Minikin, A., Lohmann, U., Ovar- 
lez, J., and Stohl, A.: Freezing thresholds and cirrus cloud forma- 
tion mechanisms inferred from in situ measurements of relative 
humidity, Atmos. Chem. Phys., 3, 1791-1806, doi: 10.5 194/acp- 
3-1791-2003, 2003. 

Herzog, A. and Vial, F.: A study of the dynamics of the equato- 
rial lower stratosphere by use of ultra-long-duration balloons, J. 
Geophys. Res., 106, 22745-22761, 2001. 

Heymsfield, A. J. and Platt, C. M.: A parameterization of the parti- 
cle size spectrum of ice clouds in terms of the ambient tempera- 
ture and ice water content, J. Atmos. Sci., 41, 846-855, 1984. 

Heymsfield. A. J. and Sabin, R. M.: Cirrus crystal nucleation by 
homogenous freezing of solution droplets., J. Atmos. Sci., 46, 


www.atmos-chem-phys. net/1 1/3757/201 1/ 


Atmos. Chem. Phys., 11, 3757-3771, 2011 



3770 


D. Barahona and A. Nenes: Dynamical states of low temperature cirrus 


2252-2264, 1989. 

Heymsfield, A. J. and Iaquinta, J.: Cirrus crystals terminal veloci- 
ties, J. Atmos. Sci., 57, 916-938, 2000. 

Hoyle, C. R., Luo, B. R, and Peter, T.: The origin of high ice crystal 
number densities in cirrus clouds, J. Atmos. Sci., 62, 2658-2579, 
2005. 

Jensen, E. J. and Pfister, L.: Transport and freeze-drying in 

the tropical tropopause layer, J. Geophys. Res., 109, D02207, 
doi:022 10.0 1029/02003 JD004022, 2004. 

Jensen, E. J., Smith, J. B., Pfister, L., Pittman, J. V., Weinstock, 
E. M., Sayres, D. S., Herman, R. L., Troy, R. F., Rosenlof, K., 
Thompson, T. L., Fridlind, A. M., Hudson, P. K., Cziczo, D. J., 
Heymsfield, A. J.. Schmitt, C., and Wilson, J. C.: Ice supersatu- 
rations exceeding 100% at the cold tropical tropopause: implica- 
tions for cirrus formation and dehydration, Atmos. Chem. Phys., 
5, 851-862, doi: 10.5 194/acp-5-85 1-2005, 2005. 

Jensen, E. J., Pfister, L., Bui, T. V., Lawson, P, Baker, B., Mo, Q., 
Baumgardner, D., Weinstock, E. M., Smith, J. B., Moyer, E. J., 
Hanisco, T. F., Sayres. D. S., Clair, J. M. St., Alexander, M. J., 
Toon, O. B., and Smith, J. A.: Formation of large (~100 pm) 
ice crystals near the tropical tropopause, Atmos. Chem. Phys., 8, 
1621-1633, doi: 10.5 194/acp-8- 162 1-2008, 2008. 

Jensen, E. J., Pfister, L., Bui, T.-P, Lawson, P., and Baumgardner, 
D.: Ice nucleation and cloud microphysical properties in tropical 
tropopause layer cirrus, Atmos. Chem. Phys., 10, 1369-1384, 
doi: 1 0.5 1 94/acp- 1 0- 1 369-20 1 0. 20 1 0. 

Karcher, B. and Burkhardt, U.: A cirrus scheme for global cir- 
culation models, Q. J. Roy. Meteorol. Soc., 134, 1439-1461; 
doi:1410.1002/Qj,1301, 2008. 

Karcher, B. and Lohmann, U.: A parameterization of cirrus cloud 
formation: homogeneous freezing of supercooled aerosols, J. 
Geophys. Res., 107, 4010, doi:4010.1029/2001JD000470, 2002. 

Karcher, B. and Haag, W.: Factors controlling upper tro- 

pospheric relative humidity, Ann. Geophys., 22, 705-715, 
doi: 10.5 194/angeo-22-705-2004, 2004. 

Karcher, B., Hendricks, J., and Lohmann, U.: Physically 

based parameterization of cirrus cloud formation for use in 
global atmospheric models, J. Geophys. Res., Ill, D01205, 
doi:012 10.01029/02005 JD0062 19, 2006. 

Kay, J. E., Baker, M„ and Hegg, D.: Physical controls on oro- 
graphic cirrus inhomogeneity, Atmos. Chem. Phys., 7, 3771- 
3781, doi: 10.5 194/acp-7-377 1-2007, 2007. 

Khvorostyanov, V. I., Morrison, H., Curry, J. A., Baumgard- 
ner, D. G., and Lawson, P: High supersaturation and 

modes of ice nucleation in thin tropopause cirrus: simula- 

tion of the 13 July 2002 cirrus regional study of tropical 
anvils and cirrus layers case, J. Geophys. Res., Ill, D02201, 
doi:02210.01029/02004JD005235, 2006. 

Kim, Y.-J., Eckermann, S. D., and Chun, H.-Y.: An overview of the 
past, present, and future of gravity-wave drag parameterization 
for numerical climate and weather prediction models, Atmos.- 
Ocean, 41,65-98, 2003. 

Koop, T., Luo, B., Tslas, A., and Peter, T.: Water activity as the de- 
terminant for homogeneous ice nucleation in aqueous solutions. 
Nature, 406, 611-614, 2000. 

Korolev, A. V. and Mazin, I. P: Supersaturation of water vapor in 
clouds, J. Atmos. Sci., 60, 2957-2974, 2003. 

Kramer, M., Schiller, C., Afchine, A., Bauer, R., Gensch, I., Man- 
gold, A., Schlicht, S., Spelten, N., Sitnikov, N., Borrmann, S., 


de Reus, M., and Spichtinger, P: Ice supersaturations and cir- 
rus cloud crystal numbers, Atmos. Chem. Phys., 9, 3505-3522, 
doi: 10.5 194/acp-9-3505-2009, 2009. 

Lawson. R. P, Pilson, B., Baker. B., Mo, Q., Jensen, E., Pfister, 

L. , and Bui, P: Aircraft measurements of microphysical prop- 
erties of subvisible cirrus in the tropical tropopause layer, At- 
mos. Chem. Phys., 8, 1609-1620, doi:10.5194/acp-8-1609-2008, 
2008. 

Liou, K.: Influence of cirrus clouds on weather and climate pro- 
cesses: a global perspective, Mon. Weather Rev., 114, 1 167 — 
1 199, 1986. 

Lohmann, U. and Feichter, J.: Global indirect aerosol effects: a re- 
view, Atmos. Chem. Phys.. 5, 715-737, doi:10.5194/acp-5-715- 

2005, 2005. 

Magee, N., Moyle, A. M., and Lamb, D.: Experimental de- 
termination of the deposition coefficient of small cirrus-like 
ice crystals near — 50 °C, Geophys. Res. Lett., 33, L17813, 
doi:17810.1 1029/12006GL026665, 2006. 

Murray, B. J.: Inhibition of ice crystallisation in highly viscous 
aqueous organic acid droplets, Atmos. Chem. Phys., 8, 5423- 
5433, doi: 10.5 194/acp-8-5423-2008, 2008. 

Murphy. D. M. and Koop, T.: Review of the vapour pressures of ice 
and supercooled water for atmospheric applications, Q. J. Roy. 
Meteorol. Soc., 131, 1539-1565, 2005. 

Murray, B. J.. Knopf, D. A., and Bertram, A. K.: The formation of 
cubic ice under conditions relevant to Earth's atmosphere. Na- 
ture, 434, 202-205, 2005. 

Murray, B. J., Wilson, T. W„ Dobbie, S„ Cui, Z„ Al-Jumur, S. M. 
R. K., Mohler, O., Schnaiter, M., Wagner, R., Benz, S., Niemand, 

M. , Saathoff, H., Ebert, V., Wagner, S., and Karcher, B.: Hetero- 
geneous nucleation of ice particles on glassy aerosol under cirrus 
conditions, Nature Geosci., 3, 233-237, 2010. 

Peter, T.: Microphysics and heterogeneous chemistry of polar 
stratospheric clouds Annu. Rev. Phys. Chem., 48, 785-822, 
1997. 

Peter, T., Marcolli. C., Spichtinger, P., Corti. T., Baker, B., and 
Koop, T.: When dry air is too humid. Science, 314, 1399-1401, 

2006. 

Pruppacher, H. R. and Klett, J. D.: Microphysics of clouds and pre- 
cipitation 2nd ed., Kluwer Academic Publishers, Boston, MA, 
954 pp„ 1997. 

Ramanathan, V. and Collins, W.: Thermodynamic regulation of 
ocean warming by cirrus clouds deduced from observations of 
1987 El Nino, Nature, 351, 27-32, 1991. 

Sato, K.: Vertical wind disturbances in the troposphere and lower 
stratosphere observed by the MU radar, J. Atmos. Sci., 47, 2803- 
2817, 1990. 

Seinfeld, J. H.: Clouds, contrails and climate. Nature, 391, 837- 
838, 1998. 

Seinfeld, J. H. and Pandis, S. N.: Atmospheric Chemistry and 
Physics, John Wiley and Sons, New York. NY, USA, 1998. 

Selkirk, H. B., Vomel, H., Valverde Canossa, J. M., Pfister, L.. 
Diaz, J. A., Fernandez, W., Amador, J.. Stolz, W., and Peng, 
G. S.: Detailed structure of the tropical upper troposphere and 
lower stratosphere as revealed by balloon sonde observations of 
water vapor, ozone, temperature, and winds during the NASA 
TCSP and TC4 campaigns, J. Geophys. Res., 115, D00J19, 
doi:10.1029/2009jd0 13209, 2010. 

Shilling, J. E., Fortin, T. J., and Tolbert, M. A.: Depositional ice nu- 


Atmos. Chem. Phys., 11, 3757-3771, 2011 


www.atmos-chem-phys. net/1 1/3757/201 1/ 



D. Barahona and A. Nenes: Dynamical states of low temperature cirrus 


3771 


cleation on crystalline organic and inorganic solids, J. Geophys. 
Res., 1 1 1 , D 1 2204; doi: 1 2 1 01 . 1 1029/12005 JD006664, 2006. 
Spichtinger, P. and Gierens, K. M.: Modelling of cirrus clouds 
- Part 2: Competition of different nucleation mechanisms, At- 
mos. Chem. Phys., 9, 2319-2334, doi: 10.5 194/acp-9-23 19-2009, 
2009a. 


Spichtinger, P. and Gierens, K. M.: Modelling of cirrus clouds - 
Part lb: Structuring cirrus clouds by dynamics, Atmos. Chem. 
Phys., 9, 707-719, doi: 10.5 194/acp-9-707-2009, 2009b. 

Wise, M. E., Baustian, K. J., and Tolbert, M. A.: Inter- 

nally mixed sulfate and organic particles as potential ice 
nuclei in the tropical tropopause region, 107, 6693-6698, 
doi:6610.1073/pnas. 0913018107, 2010. 


www.atmos-chem-phys. net/1 1/3757/201 1/ 


Atmos. Chem. Phys., 11, 3757-3771, 2011 



